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Radio lobes and X-ray hot spots in the microquasar S26 
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ABSTRACT 

We have studied the structure and energetics of the powerful microquasar/shock- 
ionized nebula S26 in NGC 7793, with particular focus on its radio and X-ray proper- 
ties. Using the Australia Telescope Compact Array, we have resolved for the first time 
the radio lobe structure and mapped the spectral index of the radio cocoon. The steep 
spectral index of the radio lobes is consistent with optically-thin synchrotron emission; 
outside the lobes, the spectral index is flatter, suggesting an additional contribution 
from free-free emission, and perhaps ongoing ejections near the core. The radio core is 
not detected, while the X-ray core has a 0.3-8 keV luminosity ~6x 10 36 erg s _1 . The 
size of the radio cocoon matches that seen in the optical emission lines and diffuse 
soft X-ray emission. The total 5.5-GHz flux of cocoon and lobes is « 2.1 mJy, which 
at the assumed distance of 3.9 Mpc corresponds to about 3 times the luminosity of 
Cas A. The total 9.0-GHz flux is w 1.6 mJy. The X-ray hot spots (combined 0.3-8 
keV luminosity w 2 x 10 37 erg s _1 ) are located « 20 pc outwards of the radio hot 
spots (i.e., downstream along the jet direction), consistent with a different physical 
origin of X-ray and radio emission (thermal-plasma and synchrotron, respectively). 
The total particle energy in the bubble is ~ 10 53 erg: from the observed radio flux, we 
estimate that only ~ a few 10 50 erg are stored in the relativistic electrons; the rest is 
in protons, nuclei and non-relativistic electrons. The X-ray-emitting component of the 



gas in the hot spots contains ~ 10 erg, and ~ 10 erg over the whole cocoon. We 
suggest that S26 provides a clue to understand how the ambient medium is heated by 
the mechanical power of a black hole near its Eddington accretion rate. 
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1 INTRODUCTION 

The basic physical model for radio lobes in FRII radio 
galaxies is based on a pair of relativistic, collimated jets 
emerging from the active black hole (BH). As the jet 
interacts with and is decelerated by the ambient (inter- 
stellar or intergalactic) medium, a reverse shock propa- 
gates inwards into the ejected plasma. After crossing the 
reverse shock, the jet material inflates a cocoon of hot 
gas, which is less dense but much overpressured with re- 
spect to the undisturbed medium. Thus, the cocoon ex- 
pands supersonically, driv ing a forward shock (bow shock) 
into t he ambient medium llScheuerll 1974c iBlandford fc Reed 
1 1974 iBegelman et al.l Il984l ; iRawlings fc Saunders! ll99lF 
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Kais er fc Alexander! Il997h . The cocoon and lobes are 
the main sources of optically-thin (steep spectrum) syn- 
chrotron radio emission, while we expect optically-thick 
(flat-spectrum) radio emission from the jet near the core. 
A radio- and sometimes X-ray-luminous hot spot is usu- 
ally found at the reverse shock, at the end of the jet. This is 
where most of the bulk kinetic energy of the jet is transferred 
to thermal ions, and to a non-thermal population of ultra- 
relativistic electrons, which cool via synchrotron and syn- 
chrotron self-Compton emission. Non-thermal X-ray emis- 
sion at the hot spot position may be due to synchrotron 
and synchrotron self-Compton emission. Optically-thin ther- 
mal plasma X-ray emission may come instead from the hot, 
shocked ambient gas between the reverse shock and the bow 
shock; in this case, the peak of the thermal X-ray emission 
will appear just in front of the radio hot spots. 
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There is a scale invariance between the jet emission pro- 
cesses in microquasars (powered by stellar-mass BHs) and in 
AGN/quasars (powered by supermassive BHs). There is also 
at least one important difference: microquasars are mostly 
located in a relatively low-pressure medium as compared to 
the medium around A GN, when s caling of the jet thrust 
is taken into account (|Hein j [20021 ) . As a consequence, we 
expect to see fewer, dimmer cocoons and radio lobes in mi- 
croquasars than in the most powerful AGN and quasars; 
however, the linear sizes of those microquasar cocoons and 
jets can be up to 1000 times larger than in radio galaxies, 
scaled to their respective BH masses. There is also evidence 
that some microquasars are located inside low-density cav- 
ities, compared with the undisturbed interstellar medium 
i|Hao fc Zhandl2009l ). 

So far, our knowledge of the interaction of microquasar 
jets with the interstellar medi um has largely relied on the 
Galactic microquasar SS433 (|Fabrikal [2004 ) and its sur- 
rounding synchrotron-emitting nebula W50 (size ~ 100 x 50 
pc). A mildly relativistic (vj = 0.27c), precessing jet acts 
as a sprinkler that inflates "ear-like" lobe structures, pro- 
truding from the more spherical W50 nebula. Most of 
the jet power [~ 10 39 erg s _1 ) is dissipated in the lobes 
l|Begelman et al.l Il980l ). Faint evidence of the interaction 
of relativistic jets with the interstellar medium has been 
found in a few other, le ss powerful Galac tic microquasars: 
for example Qyg X -l (|Gallo et al.l l2005h, GRS 1915+10 5 
jKaiser et al1|2004l). XTEJ1550 -564 (|Corbel et~aH l2002t) . 



H1743-322 JCorbel et alj|2005l), and around the neu tron 



star Sco X-l ( Fomalont. Geldzahler fc B radshaw 200l|). On 
a larger scale, huge (size > 100 pc) ionized nebulae have 
been found around several ultraluminous X-ray sources 
(ULXs) in nearby galaxies (fPakull fc Mirionil |2002|. | 2003 |: 



Roberts et al. 2003; Pakull ct al. 2006; Pakull & Grisc 2008; 



Grise et al.ll2008l : iFeng fc Kaaretl 12008). Such nebulae emit 



optical lines typical of shock-i onized gas, and in a few cases, 
synchrotron radio emi ssion ([Miller et all 120051 : [§ oria et al.l 
l2006l : lLang et al.ll2007l ). The derived ages (> a few 10 5 yr) 
and energy content (~ 10 52 -10 53 erg) are too large for or- 
dinary supernova remnants, and suggest jet/wind inflation 
with a mechanical power ~ 10 39 -10 40 erg s 1 , comparable 
with the X-ray luminosities l| Pakull et al J 1200a ). However, 
no direct X-ray or radio evidence of a collimated jet has been 
found in ULX bubbles so far. On the other hand, X-ray lu- 
minous sources may be only a s ubset of non-nuclear B Hs at 
very high mass accretion rates. IPakull fc Grisi | 2008) pro- 
posed that ionized bubbles might also be found associated 
with BHs that appear X-ray faint, either because their ra- 
diative emission is collimated away from our line of sight, 
or because they are transients and currently in a low/off 
accretion state, or because they channel most of their accre- 
tion power into a jet even at near-Eddington mass accretion 
rates. 



2 THE MICROQUASAR S26 IN NGC 7793 

A spect acular example of such syst ems was recently dis- 
covered jPakull. Soria fc MotchhoiOl . henceforth PSM10) in 
the outskirts of the Sculpt or galaxy NGC 7793 (Figu re 1), 
at a distance of 3.9 Mpc ([Karachentsev et al.l 12003). The 
radio/optical nebula S26 was originally classified as a su- 



Figure 1. Location of S26 inside its host galaxy NGC 7793, from 
public-archive BVR images taken from the Cerro Tololo Inter- 
American Observatory (CTIO) 1.5-m telescope on 2001 October 
18. S26 is inside the circle (radius of 15") in the north-eastern 
corner of the galaxy. 



pernova remnant candidate CBl air fc Longl [T 997i ) ; the high 
[S II] A6716, 6732/Ha flux ratio indicates the presence of 
shock-ionized gas. The optical radial velocity of S26 agrees 
with that of NGC 7793, ruling out a chance superposition 
of a background AGN. A radio spectral index consistent 
with optically -t hin s ynchrotron emission was reported by 
iPannuti et alj (|2002l ). and the emitting region appeared 
clearly extended and elongated. However, the spatial resolu- 
tion was too low to reveal details of its internal structure. A 
faint X-ray source was di scovered to be associa ted with S26 
in ROSAT observations jRead fc Pietschlll999l), but it was 



Pakull fc Gris. 



i i[2008l ) 



dis- 



unresolved. Using Chandra data, 
covered that the X-ray emission is resolved into three sources 
that are perfectly aligned and match the extent of the ma- 
jor axis of the radio and optical nebulae. Those sources have 
been interpreted as the core (at the X-ray binary position) 
and the X-ray hot spots (where the jet interacts with the 
ambient medium). 

From optical spectroscopic observations, PSM10 deter- 
mined the expansion velocity, density and temperature of 
the line-emitting gas in the bubble, and discovered that the 
mechanical power of the central BH is ~ a few x 10 40 erg s _1 ; 
this suggests accretion rates similar to those required for the 
most luminous ULXs. PSM10 showed that the jet power is 
orders of magnitude higher than both the X-ray luminosity 
and the value one would derive from the radio luminosity; 
they argued that most of the jet power is transferred to 
non-relativistic protons and nuclei rather than non-thermal 
relativistic electrons. 

In this paper, we present the initial results of our radio 
study, showing for the first time the resolved lobe struc- 
ture and measuring the spectral index variations across the 
source. We discuss the origin of the radio emission and the 
implied jet power. We compare the radio, X-ray and op- 
tical maps of the nebula, determining the positions of the 
radio and X-ray hot spots and of the core, and we provide 
a more detailed spectral analysis and interpretation of the 
X-ray properties. We then summarize the energy budget of 
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Figure 2. Top panel: ATCA map at 5.48 GHz, with in- 
tensity scale in Jy beam -1 . Contour intensity levels are 
(0.40,0.84,1.29,1.73,2.17,2.61,3.06,3.50) xlO -4 Jy beam" 1 . 
The rms noise level is 0.085 X 10 — 4 Jy beam -1 . Overplot- 
ted green circles mark the positions of the X-ray core and 
hot spots (0"6 radius for the core, 0'.'8 for the hot spots), 
from Chandra. Bottom panel: ATCA map at 9.02 GHz, with 
intensity scale in Jy beam -1 . Contour intensity levels are 
(0.39,0.65,0.91,1.17,1.44,1.70) xlO -4 Jy beam -1 . The rms 
noise level is 0.135 x 10 -4 Jy beam -1 . 
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Figure 3. Map of the radio spectral index, inferred from the ratio 
of the 5.48 GHz and 9.02 GHz maps, where the latter has been 
tapered to the resolution of the former. Green circles mark the 
positions of the X-ray core and X-ray hot spots. 

this system, quantifying the fraction of energy stored in rel- 
ativistic electrons and in the X-ray emitting gas. 



3 OBSERVATIONS 
3.1 Radio observations 

We observed S26 on 2009 August 6 and 7 with the Aus- 
tralia Telescope Compact Array (ATCA). Simultaneous 5.5 
and 9 GHz observations were carried out with the Compact 
Array Broadband Backend (CABB); the bandwidth at each 
frequency is about 2 GHz. The array configuration was 6D, 
with minimum and maximum baselines of 77 and 5878 m, 
respectively. The total integration time on-source was 13.3 
hours; the data for antenna 6 at 9 GHz were lost during part 
of the second observing session due to technical problems. 
B1934— 638 was used as the primary calibrator, while our 
secondary calibrator was B2357— 318. 

We reduced and imaged the data with MIRIAD 
jSault et al.lll995l ). After flagging bad data, the effective 
frequencies of the two bands are 5.48 and 9.02 GHz. We 
tried different values of Briggs' robust weighting parame- 
ter |Briggdll1)9"5T ); we found that a value of 0.0 provides 
a good balance between sidelobe suppression and sensitiv- 
ity at both frequencies. Because of the wide bandwidths, 
we used the multi-frequency dec onvolution algorithm MF- 
CLEAN l|Sault fc Wieringal Il994l ) . The CLEANed, primary- 
beam-corrected images are shown in Figure 2; the angu- 
lar resolutions are 3" 54 x l'.'38 (position angle 1.1°) and 
2'.'67 x 1.08 (position angle 6.8°) at 5.48 and 9.02 GHz, re- 
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GHz data. We created a two-point spectral index map (Fig- 
ure 3), where the sign of the index is defined such that the 
specific flux Si, ~ v a . For the tapered 9.02 GHz data, we 
found that a robust weighting parameter of 0.5 provides the 
best compromise between residual sidelobe contamination 
and sensitivity to the low-surface-brightness extended emis- 
sion that is clearly visible in the 5.5 GHz map. 



Figure 4. Top panel: Chandra/ ACIS color map of S26 (smoothed 
with a 1" Gaussian core), with ATCA 5.48-GHz intensity con- 
tours. The colour coding is: red: 0.2-1 keV; green: 1-2 keV; blue: 
2-8 keV. Contour intensity levels are as in Figure 2. Bottom panel: 
Chandra/ ACIS X-ray colour map (same colour coding), with 
ATCA 9.02-GHz radio intensity contours superimposed. Contour 
intensity levels are as in Figure 2. Cf. Fig. 1 in PSM10, which 
overplots Ha emission contours over the same Chandra image. 



spectively. In the vicinity of the microquasar, the rms noise 
levels are 8.5/xJy beam -1 (5.48 GHz) and 13.5/xJy beam -1 
(9.02 GHz). We estimate that the internal calibration un- 
certainty is ~ 2 per cent at both frequencies. 

We also tapered the 9.02 GHz data so that the reso- 
lution and beam position angle matched those of the 5.48 



3.2 X-ray observations 

NGC7793 was observed with Chandra/ ACIS-S3 on 2003 
September 6 (Obs ID 3954). The live time was 48.9 ks. We 
retrieved the data from the public archives (processed with 
ASCDVER=7.6.8), and analysed them with standard imag- 
ing and spectroscopic tools such a s psextract in the dat a 
analysis system CIAO Version 4.0 (jFruscione et al.l l2006t ) . 
We modelled the X-ray spectra with XSPEC Version 12.0 
|Arnaudlll996h . Luckily, the roll angle of the Chandra ob- 
servation was such that S26 was located rather close to the 
S3 aimpoint (less than 1' away), giving us a narrower point 
spread function. 



4 MAIN RESULTS 
4.1 Radio results 

The most important new result of our ATCA study is 
that we have resolved the spatial structure of the radio- 
emitting nebula. Most of the emission comes from two ra- 
dio hot spots and surrounding lobes, with a fainter but 
clearly identified cocoon encompassing the m (Figure 2, top 
pane l) . This is the textbook structure (e.g.. lBegelman et al.l 
1984) of FR II- type powerful radio galaxie s (e.g., Cygnus A: 
Carilli fc Barthelll 19961 : IWilson et aUl2006t ). The radio struc- 
ture is aligned with the jet axis suggested by the three X- 
ray sources, confirming this interpretation. The position of 
the northern radio hot spot is RA = 23 h 57 m 59 s .58, Dec 
= -32°33'13'.'6 (with an uncertainty of w 0'.'2). The posi- 
tion of the southern hot spot is RA = 23 h 58 m 00M5, Dec 
= — 32°33'25'.'0. Thus, the projected distance between the 
radio hot spots is (13'.'5 ± 0'.'3) » 250 pc. We interpret the 
radio hot spots as the reverse shocks (Mach disks) at the 
ends of the jets. 

At 5.5 GHz, the peak intensity in the southern lobe 
is ~ 0.37 mjy beam -1 ; in the northern lobe, « 0.21 mjy 
beam -1 ; the total flux in the lobes and cocoon is « 2.1 mjy 
(Table 1), that is » 3 times the luminosity of Cas A. From 
the untapered map at 9 GHz, we obtain a peak intensity in 
the southern lobe ~ 0.19 mjy beam -1 ; in the northern lobe, 
« 0.11 mjy beam -1 . The total flux at 9 GHz is « 1.6 mjy 
(Table 1). The spectral index in the lobes is, on average, 
~ —0.7 to —0.6; it appears to be flatter (~ —0.4 to 0) across 
most of the cocoon, and inverted (« to 0.4) at the base 
of the jets, on either side of the X-ray/optical core (Figure 
3). We estimate a la uncertainty for a of « 0.12 near the 
southern radio hot spot, sa 0.19 near the northern radio hot 
spot, and « 0.5-0.6 in the rest of the cocoon, where the 
emission is much fainter. Thus, the existence of a complex 
spatial structure for the spectral index is at this stage still 
an intriguing speculation that has to be tested with deeper 
observations. 
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Figure 5. Left panel: unsmoothed Chandra color map, with the same colour coding as in Figure 4. Middle panel: greyscale, continuum- 
subtracted Ha emission, with ATCA 5.5-GHz intensity contours. Red circles mark the positions of the X-ray core and hot spots (0'.'6 
radius for the core, 0'.'8 for the hot spots). The public-archive Ha image was taken from the CTIO 1.5-m telescope on 2001 October 
18. Right panel: greyscale, continuum-subtracted Hell A4686 emission, with ATCA 9.0-GHz intensity contours and X-ray core/hot spot 
positions. The image was taken with the FORS1 camera on the ESO Very Large Telescope on 2002 November 1; we smoothed it with a 
0'.'6 Gaussian core. See PSM10 for a discussion of the bright optical core. 



Structure 


Sb.5 (mJy) 


S 9 .0 (mJy) 


Total 


2.1 ±0.1 


1.6 ±0.1 


S lobe 


0.95 ± 0.05 


0.80 ±0.05 


N lobe 


0.45 ± 0.03 


0.28 ±0.03 



Table 1. Observed radio fluxes for the whole nebula, and for the 
two lobes, from our ATCA observations. The 9.0-GHz fluxes are 
measured from the tapered map. 



4.2 X-ray results 

The key feature of this system is the aligned triplet of point- 
like sources (Figur es 4, 5), which we int erpret as the X-ray 
core and hot spots (|Pakull fc G rise 200<|). The X-ray core is 
located at RA = 23 ,l 57 ra 59 s .94, Dec = -32°33'20'.'9 (with 
an uncertainty of « 0?2). It has a hard spectrum (power- 
law photon index r = 1.4 ± 0.6), consistent with a BH in 
the low/hard state (Rcmillar d fc McClintockll2006l ). and an 
emitted luminosity Lo.3-8 ~ 6 x 10 36 erg s _1 (Table 2 and 
Figure 6). The X-ray core coincides, within the astromet- 
ric uncertainties, with a point-like optical source with Hell 
A4686 emission (Figure 5). The 90% uncertainty circle of 
the ACIS-S3 absolute position has a radius of 07 4j, and the 
uncertainty of the optical images is « 0'.'3. 

The two hot spots have a much softer spectrum (Fig- 
ures 4, 6), and are well fitted (Table 3) by a 2-component 
ray m ond- smith thermal plasma model l|Ravmond fc Smithl 
Il977l ) with feTi w 0.3 keV and kT 2 » 0.9 keV, and neg- 
ligible intrinsic absorption (Cash statistics = 10.4 over 13 
dof for solar abundances; = 9.6 over 13 dof for 1/4-solar 
abundances). The emitted X-ray luminosities are Lo.3-8 ~ 



5 x 10 36 erg s _1 and L .-. 



11 x 10 36 erg s" 1 for the 



northern and southern hot spot, respectively (similar to the 



http : / /cxc. harvard .edu /cal/ASPECT /celmon / 



ratio of radio luminosities). In general, a choice of low metal 
abundances give better fits than solar abundances, but the 
signal-to-noise is not high enough to contrain this param- 
eter. Other, more complex thermal plasma models such as 
mekal, vmekal, equil and nei also give similar sets of best- 
fitting parameters; they all require at least two tempera- 
ture components. However, the sedov thermal plasma model 
(jBorkowski. Lverlv fc Revnolds!l200ll ) gives a good fit (Cash 
statistics = 11.9 over 15 dof for 1/4-solar abundances) with 
only one temperature component, at kT « 0.52 keV (Ta- 
ble 4); the ionization age (r parameter) of the sedov model 
is consistent with the characteristic age multiplied by elec- 
tron density in the S26 bubble. Given the small number of 
counts in the Chandra spectrum, none of the more com- 
plex thermal-plasma models can provide an improvement 
over the simpler raymond- smith model. Simple or broken 
power-law models do not give acceptable fits (Cash statis- 
tics = 25.5 over 15 dof); moreover, they would require an 
unphysically steep slope (F « 6) combined with high in- 
trinsic column densities (Nh ~ 5 x 10 21 cm -2 ). We con- 
clude that the hot spot spectra are not dominated by syn- 
chrotron or synchrotron self-Compton emission. We inter- 
pret them as optically-thin thermal plasma emission from 
hot, shocked gas, probably located between the reverse and 
forward shocks. The radius of the X-ray hot spots is < 1". 
From the combined volume of the two hot spots and their 
emission measures (Table 3), we estimate a hot gas density 
> 1 cm and a mass ~ 10 36 g (see also PSM10). 



We also find faint X-ray emission projected over the 
surface of the cocoon (Figure 5), with slightly softer colours 
than the hot spots; however, the number of detected counts 
is too low for detailed temperature comparisons. A single- 
temperature bremsstrahlung fit suggests kT = 0.5 ±0.1 keV 
(Table 5 and Figure 7). It is also impossible to determine at 
this stage whether the X-ray emitting gas is filling the whole 
cocoon, or is confined to an outer shell. If the hot gas fills a 
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Parameter Value 



Aff.Qal 

r 

K 


1.2 X IO 20 cm" 2 (fixed) 
(4.3t 2 ;J) x 10 J7 


/0.3-s 
£0.3-8 


(3.4±g |) X 10 -15 erg cm" 2 s" 1 
(6.2ti'g) x 10 36 erg s" 1 


C-statistic 


4.56 (6 dof) 



Table 2. Best-fitting spectral parameters for the X-ray emission 
from the core. The XSPEC model is wabsQ^i *powerlaw. Adding 
intrinsic absorption does not improve the fit. Errors are 90% con- 
fidence level for 1 interesting parameter. Here, and in the follow- 
ing tables, the line-of-sight column density to NGC 7793 is from 
iKalberla etafl J2005| V 

significant fraction of the cocoon, its characteristic density 
is ~ 0.1 cm' 3 and the total mass ~ 10 37 g. 

The position of the northern X-ray spot is RA = 
23 h 57 m 59".56, Dec = -32°33'12'.'6 (with an uncertainty of 
« Of! 2). The position of the southern X-ray spot is RA = 
23 ,l 58 m 00 a .20, Dec = -32°33'25'.'7. Therefore, the projected 
distance between the X-ray hot spots is (15'.' 4 ± 0'.'3) ~ 290 
pc, slightly larger than the projected distance between the 
radio hot spots; that is, each X-ray hot spot appears pro- 
jected m 15-20 pc downstream from the corresponding radio 
hot spot. In fact, the X-ray and radio hot spot positions in 
the southern lobe coincide with two distinct brightness peaks 
in the optical images, particularly Ha and V band (Figures 
5, 8). This provides another argument in favor of two sepa- 
rate emission mechanisms for the radio and X-ray hot spots. 
The core is along the line of the hot spots, but not symmet- 
rically located between them. It is projected w 6" from the 
southern spot (which is about twice as luminous in every 
band) and ~ 9'.' 5 from the northern spot. This may perhaps 
be due to a higher density of the interstellar medium south 
of the core. 



Parameter 


Value 


z 

fcTi 
Ni 
kT 2 
N 2 


1.2 X 10 20 cm" 2 (fixed) 
< 1.0 X 10 21 cm" 2 
1 (fixed) 
°.26lS:i£ keV 
(2.3+ 2 ; 9 6 ) x 10-6 
0.96±°; 3 2, keV 
(l.9±g"J) x 10- 6 


/0.3-8 
£0.3-8 

EM(0.26 keV) 
EM(0.96 keV) 


( 8 - 8 tos>) x 10-15 er § cm ~ 2 s_1 
(l.7lo J) x 1037 er S s_1 
) X 10 59 cm" 3 
(3.5ti; 3 ) x 10 59 cm" 3 


C-statistic 


10.35 (13 dof) 


Wff.Gal 

N H 
Z 

fcTi 
Ni 
kT 2 
N 2 


1.2 x 10 20 cm" 2 (fixed) 
< 1.1 X 10 21 cm" 2 
0.25 (fixed) 
0.29±°;°8 keV 
(7.3± 3 2 i 2 ) x IO" 6 
0.90+q 1 ; 32 keV 
(4.9+ 9 ; 2 ) x IO" 6 


/o.3 — 8 
^0.3-8 

EM(0.29 keV) 
EM(0.90 keV) 


(9.1+S o) X 10 -15 erg cm -2 s" 1 

I (Jul '-' 

(l.8±o |) X 10 37 erg s" 1 
(l3.3t 22 6 2 ) x 10 59 cm- 3 
(8.9+J 6 5 7 ) x 10 59 cm" 3 


C-statistic 


9.55 (13 dof) 



Table 3. Best-fitting spectral parameters for the (combined) hot 
spot X-ray emission, assuming solar and 1/4-solar abundances. 
The XSPEC model is wabsQ !i i*wabs*(ray+ray). Errors are 90% 
confidence level for 1 interesting parameter. 



5 DISCUSSION 

5.1 Energetics of the bubble 

We have presented radio and X-ray results from our multi- 
band study of a powerful non-nuclear BH in NGC 7793, and 
of its surrounding shock-ionized cocoon (see PSM10 for a 
discussion of the evidence for shock ionization from the op- 
tical emission lines). The system was ori ginally classified as 
a supernova re mnant (jBlair fe Londl 19971 ). In that scenario, 
Asvarov (2006) showed that an input energy « 5 x 10 52 erg 
was required to explain its size and radio luminosity, well be- 
yond the energy that can be supplied by an individual super- 
nova. Based on the clear radio and X-ray evidence for a colli- 
mated jet pair (lobes, hot spots), PSM10 showed that such a 
large amount of energy has been supplied by the BH over the 
lifetime of the bubble (characteristic age « 2 x 10 5 yr). The 
core is seen as a faint point-like X-ray source, consistent with 
a stellar-mass BH in the low/hard state, and a point-like op- 
tical source, consistent with an OB donor star (possibly a 



Wolf-Rayet: PSM10). It is undetected in the radio bands, to 
a 3<t upper limit ~ 0.03 mjy. This is unsurprising: if the BH 
lies in the fundamental plane l|Merloni. Heinz fc Di Matted 
l2003l ; lKording. Falcke fc Corbelll 20061) with a mass ~ 10M© , 
we expect a core radio flux ~ 0.01/xJy (yL v ~ 10 30 erg 
s _1 ), like from a common-or-garden low/hard state micro- 
quasar at a distance of 3.9 Mpc. On the other hand, the 
impact of this BH onto the surrounding interstellar medium 
is all but common. The simplest explanation is that the 
core is currently in a low/hard state, three or four orders 
of magnitude fainter than its long-term average power. In 
that same canonical state, the s teady jet power V oc Ly 5 
iFender, Belloni fc Gallol 12004 llvlalzac. Merloni fc Fabianl 



|2004 iFender. Gallo fc Jonkeij 120031 '). The normalization of 
this relation has an uncertainty of almost two orders of 
magnitude, but is constrained enough to suggest 10 36 erg 
s_1 ^5 10 38 erg s _1 , also much lower than the inferred 
long-term average. However, there may be alternative expla- 
nations for the apparent faintness of the core. Perhaps the 
observed X-ray luminosity is a severe underestimate of the 
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Parameter 



Value 



Parameter 



Value 



Aff.Qal 

N H 
Z 

kr 

N 


1.2 x 10 2 ° cm -2 (fixed) 
< 2.5 X 10 21 cm- 2 
0.25 (fixed) 
0.78±° °l keV 
(SAtfl) x 10~ 6 


C-statistic 


28.9 (15 dof) 


Nrr 

z 

kT a 
kT h 

T 

N scd 


1.2 X 10 20 cm" 2 (fixed) 
(fixed) 
0.25 (fixed) 
0.52+° \l keV 
kT a (fixed) 

(l-9lj V) x 1Ql2 s cm ~ 3 
f8.7 +2 S) x 10~ 6 

V — 2.8/ 


/o.3-8 
^0.3-8 

EM 


(9-Oli g) X 10" 15 erg cm" 2 s" 1 

(l.slg x 1037 er e s_1 

(l5.9lg i) x 10 59 cm" 3 


C-statistic 


11.88 (15 dof) 



Table 4. Alternative spectral models for the combined hot spot 
X-ray emission. A single-temperature wabsQ a i*wabs*ray model 
does not produce acceptable fits. However, a single-temperature 
wabsQ^i *wabs*sedov model results in a fit as good as those with 
two-temperature ray mond- smith models (Table 3). Errors are 
90% confidence level for 1 interesting parameter. 



true X-ray luminosity, if most of the direct emission is ab- 
sorbed and/or beamed awa y from our line of sight (as i t has 
been suggested for SS 433: IMedvedev fc Fabrikall201Gi ) and 
we are only seeing a scattered component. Or perhaps the 
system is not in the canonical low/hard state, but in some 
other unclassified state with V 3> Lx ■ 

The faint core is in stark contrast with the large, bright 
nebula, visible in all bands with a similar size and shape 
(Figure 5), and a conservatively estimated volume « 10 62 
cm -3 , assuming a prolate spheroid with a major axis « 280 
pc and minor axes « 130 pc (based on the projected distance 
between the hot spots along the major axis, and the width 
of the radio nebula at 5.5 GHz in the transverse direction). 
In fact, Ha images (PSM10) may suggest an even larger 
size, « 340 x 170 pc. Thus, the volume-averaged shell radius 
R s ~ 100 pc. Its characteristic size is an order of magnitude 
larger than the jet driven bubble around Cy g X-l, which 



has a n estimated jet power ~ 10 37 erg s 1 l|Russell et al.l 



12007ft . It is a factor of two larger (allowing for distance un- 
certainties) than the radius of the SS4 33/W50 nebula, with 
an estimated jet power ~ 10 39 erg s" 1 (IMedvedev fc Fabrikal 
l20ld : lFabrika|[200l : iMarshall. Canizares fc Schulj|2002l ). ~ 
PSM10 determined a mechanical power V ~ 5 x 10 40 
erg s _1 for S26, using the well-known self-similar solution to 
the conservation of mass, momentum a nd energy equations 
(Equations 17-22 in l Weaver et al.lll977f ). in which the radius 
of the swept-up shell 7? s « 0.76 V 1/5 t 3/5 pg 1/5 . They mea- 
sured the expansion velocity from the half-width at zero- 



1 * rl ,tjal 


1.2 X 10 20 cm -2 


(fixed) 


z 


0.25 


(fixed) 


feTi 


0.29 keV 


(fixed) 


Ni 


(3-5±?:S) 


x 10~ 6 


kT 2 


0.90 keV 


(fixed) 


N 2 


(2-i±i;l) 


x 10~ 6 



jb.3-8 


( 4 - 2 to<i) x 10 ~ 15 er § cm ~ 2 s ~ 


-1 


^0.3-8 


(7.§±\\) X !0 36 erg s" 


-l 


EM(0.26 keV) 


(6A+1 1) X 10 59 on" 


-3 


EM(0.96 keV) 


(3.81 2 ; 9 ) x 10 59 cm" 


-3 



C-statistic 


11.28 (17 dof) 


kT br 
Kbx 


1.2 X 10 20 cm" 2 (fixed) 
OATt^l keV 
(7.5^; 9 ) x 10-6 


/0.3-8 
Lo.3-8 
EM 


(5-llo I) x lO" 15 erg cm' 2 s" 1 
(9.3ti;g) x 10 36 erg s" 1 
(4-5+ 3 fj X 10 60 cm- 3 


C-statistic 


9.77 (17 dof) 



Table 5. Best-fitting spectral parameters for the X-ray emission 
from the cocoon (not including the hot spots). The XSPEC mod- 
els are wabsQ^i *ray and wabsQ & \ *brems. Adding intrinsic absorp- 
tion does not improve the fit. Errors are 90% confidence level for 
1 interesting parameter. 



intensity of the optical emission lines (v w 250 km s _1 ) 
and from shock-ionization models of the Hell A4686/H/3 
flux ratio (v ~ 275 km s _1 ). This implies a characteristic 
age t = (3J? s )/(5w C xp) ~ 2 x 10 yr. The hydrogen num- 
ber density of the interstellar medium into which the bub- 
ble expands was estimated as no w 0.7 cm -3 (PSM10), 
from a comparison of the observed Ha emission with the 
intensity of a fully rad iative shock dDopita fc Sutherland! 
1 19961 : | Pakull et alj|20(i3 ); this corresponds to a mass den- 
sity po ~ fJ.m p no ~ 1.6 x 10~ 24 g cm -3 (taking the mean 
atomic weight /i = 1.38). The swept-up mass in the expand- 
ing shell is (47r/3)pi? 3 ~ 2 x 10 38 g, carrying a kinetic energy 



~ (15/77)Pi ~ 6 x 10 j2 erg. The energy content of the ther- 
mal gas between the reverse shock and the swept-up shell is 
E = (5/ll)Vt w 10 53 erg. 

The only shock-ionized nebulae of comparable size 
and energy content in the local universe are those around 
ULXs such as Holmberg IX X-l, NGC 1313 X-2, IC 342 X-l 
llPakull fc Mirionill2002t|Pakull fc Grisell2008l : [Ramsey et all 
l2006l : lFeng fc Kaaretll2008l ). One difference is that, unlike all 
previously known ULX bubbles, S26 has clear evidence of 
collimated jets. We do not know the relative distribution of 
mechanical power between the collimated jets and perhaps a 
more spherically symmetric wind (for example an accretion 
disk wind); however, the elongated structure of the nebula 
and the presence of bright lobes and hot spots suggests that 
the jet carries most of the power. Another difference is that 
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Figure 6. Chandra/ ACIS spectra of the combined hot spot emis- 
sion (red datapoints) and core emission (blue datapoints), fitted 
with a single-temperature sedov thermal-plasma model and with 
a power-law model, respectively. See Tables 2, 3, 4 for the best- 
fitting parameters. The sedov model fit illustrated here is statis- 
tically equivalent to the two-temperature raymond- smith model 
fit illustrated in Fig. 2 of PSM10. 




Figure 7. Chandra/ ACIS spectrum of the cocoon emission, fitted 
with a Ri 0.5 keV bremsstrahlung model. See Table 5 for the best- 
fitting parameters. 



in ULX bubbles, the central BH is (by definition) X-ray lu- 
minous, with an apparent X-ray luminosity ~ 10 40 erg s _1 , 
similar to the mechanical power. S26 may be an example 
of ULX bubble where the central BH is currently in a low 
and/or jet-dominated state. In that respect, S26 is analogous 
to (but two orders of magnitude more energetic than) the 
shock-ionized bubble around the very massi ve but only mod- 
erate ly luminous non-nuclear BH IC 10 X-l l|Prestwich et al.l 
l2007ft . 



5.2 Synchrotron and thermal plasma emission 

The detection of a radio cocoon with bright radio hot spots 
shows that some of the input mechanical power goes into 
synchrotron-emitting, relativistic electrons. As expected, the 
radio spectrum is steep in the lobes and flat or inverted near 



the position of the core (Figure 3); this may be evidence of 
recent or continuing ejection activity (although the core it- 
self is undetected in the radio). In addition, there is radio 
emission from the cocoon region outside the lobes, with a 
specific flux ~ 0.7 mjy at 5.5 GHz. There is circumstatial 
evidence that this emission has a rather flat spectrum, cer- 
tainly flatter than in the lobes (Figure 3). This is difficult 
to reconcile with a scenario where the synchrotron-emitting 
electrons in the whole cocoon are backflowing from the lobes. 
In that case, the spectral index in the rest of the cocoon 
would be even steeper — as we see for example in Cygnus A 
l|Carilli fc Barthellll996l ). We suggest that the extended ra- 
dio emission in the cocoon outside the lobes may have a sig- 
nificant contribution from (flat-spectrum) free-free emission, 
from the same thermal gas responsible for the optical recom- 
bination lines. For the characteristic temperature ~3x 10 4 
K estimated by PSM10, the ratio of the H/3 and free-free 
radio emissivity is jga /.75.5GHz ~ 2 x 10~ 10 erg cm" 2 s -1 
Jy -1 (|Caplan fc Deharven3ll98rj . their Appendix A). The 
H/3 luminosity of S26 is « 10 38 erg s" 1 (PSM10): thus, we 
expect a free- free radio flux w 0.3 mjy at 5.5 GHz. This 
is negligible in the lobes, compared with the synchrotron 
component, but may be significant in the region outside the 
lobes, and may explain the rather flat spectral index there. 

What fraction of the mechanical power is transferred to 
relativistic electrons? If we combine the self-similar model of 
cocoon expansion with the standard synchrotron emissivity, 
in the minimum-energy approximation, assuming a spectral 
index a = —0.7, we obtain (Appendix A): 



S v ~ 82(l + fc) 



-1 1.85 -r, 1.34.0.32 0.51 ,-2 — 0.7 t 

r\ P39 t B n-L dj v h mjy, 



(1) 



where r\ is the fraction of the total energy density contained 
in all relativistic species (electrons, protons and nuclei) plus 
magnetic field, (1 + fc) -1 is the fraction of relativistic particle 
energy carried by the synchrotron-emitting electrons alone, 
V39 is the jet power in units of 10 39 erg s _1 , £5 is the source 
age in units of 10 5 yr, ni is the interstellar number density in 
cm -3 , d is the source distance in Mpc, and V5 the observed 
frequency in units of 5 GHz. A specific flux a few times 
higher is expected if we assume a — —0.5 (Appendix A). 

If we assume that all the jet power is transferred to the 
relativistic electrons (that is, if we put k — and r) = 1), 
Equation (1) grossly overestimates the radio emission, for 
the measured jet power and distance of S26. (Or, conversely, 
the observed radio flux would lead us to underestimate the 
jet power if we did not know it independently). This tells 



us that (1 + fc) x r) 1 



10" 



We cannot separately 



determine fc and 77 from this simple model, but for plausible 
values of k ~ 10-100 found in cosmic rays, we estimate that 
the fraction (1 + fc) -1 x r\ of the total injected mechanical 
power carried by the relativistic electrons is ~ a few 10 -3 . 
The rest of the energy is given to protons, nuclei and non- 
relativistic electrons, and is used for heating and inflating 
the bubble, accelerating the shell of swept-up interstellar 
medium to the expansion speed ~ 250 km s~ . 

The X-ray emission from the hot spots and cocoon pro- 
vides another clue to understand the energy budget. In radio 
galaxies, X-ray hot spots are usually interpreted either as 
direct synchrotron, or synchrotron self-Compton em ission 
l|Harris fc Krawczvnskil 12002, ; lHardcastle et all 120041 ) ■ from 
the same population of electrons responsible for the radio 
hot spots, which are accelerated at the reverse shock. For 
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Figure 8. Close-up view of the southern lobe: true-colour image 
in the B, V, R bands, taken by Dr. Jifeng Liu with the Baade 
Magellan telescope on 2009 August 28 (exposure time: 200 s per 
filter). The ATCA 9.02-GHz intensity contours are overplotted 
in green; the red circles mark the position of the X-ray core and 
southern hot spot. The point-like optical core has a brightness 
B rs 23 mag, M B « -5 mag (PSM10). 



S26, we estimate from Chandra a specific flux w 10" 1 erg 
cm -2 s~ 2 keV" 1 at 1 keV from both hot spots, correspond- 
ing to Ri 4 x 10~ 32 erg cm -2 s~ 2 Hz -1 . The combined ra- 
dio emission from the radio hot spots is ~ 1 mjy ~ 10~ 26 
erg cm" 2 s" 2 Hz" 1 at 5.5 GHz. If we extrapolate the ra- 
dio flux as a straight power law, S u ~ ^"°' 7 , we would also 
expect r;4x 10" 32 erg cm" 2 s" 2 Hz" 1 at 1 keV. This of 
course requires a continuous acceleration of the most ener- 
getic electrons, so that there is no spectral break from the 
radio to the X-ray bands. So, in principle, the X-ray to ra- 
dio flux ratio is consistent with a simple synchrotron compo- 
nent. However, the X-ray spectrum tells a different story. Its 
shape and slope are not consistent with either synchrotron 
or inverse-Compton power-law models, even accounting for 
the low number of counts. We showed (Section 4 and Figure 
6) that the X-ray emission from the hot spots, with its peak 
at ~ 0.6-0.9 keV and its sharp drop above ss 1 keV, is most 
likely due to hot thermal plasma with a range of tempera- 
tures up to ~ 0.9 keV. We also showed that the X-ray hot 
spots are located m 1" further away from the core than the 
radio hot spots. This is a second argument in support of our 
claim that X-ray and radio hot spots are due to different 
physical processes. 

It is still not clear what is heating at least part of the 
X-ray emitting gas to such high temperatures, particularly 
at the hot spots. We suggest two alternative scenarios. The 
first scenario is that the X-ray emitting gas is the shocked 
interstellar medium, heated by the bow shock (advancing at 
a speed Vbe) to a temperature kT = (3/ 16) fj,m p v^, B . In this 
case, the X-ray hot spots mark the position of the bow shock 



and the radio hot spots that of the reverse shock into the 
ejecta. To produce temperatures ~ 0.9 keV (as required by 
our fits with raymond- smith and most other thermal plasma 
models in XSPEC), the bow shock velocity (that is, the ex- 
pansion velocity along the major axis) would have to be 
~ 900 km s" 1 , almost 4 times higher than the expansion 
velocity measured by PSM10 from the width of the optical 
lines. However, their slit position was almost parallel to the 
minor axis and did not include the hot spots; moreover, the 
viewing angle of the major axis is still unknown. Thus, we 
still do not know at what speed the jet heads are advancing 
into the interstellar medium, Besides, we noted in Section 
4.1 that a sedov model gives a good fit of the hot spot emis- 
sion with thermal plasma temperatures as low as « 0.5 keV, 
requiring a more plausible shock velocity ~ 300 km s" 1 . To 
sum up, we cannot yet rule out the fast bow shock scenario 
for S26. An example of an X-ray-emitting bow shock located 
ahead of the radio-emitti ng lobes can be s een in the near- 
est radio galaxy, Cen A (|Kraft et alj|2007h . An alternative 
scenario, considered more likely by PSM10, is that the bow 
shock is not advancing fast enough to produce the X-ray 
emitting gas, and the shocked interstellar gas between the 
bow shock and the contact discontinuity has already cooled 
and collapsed to a thin, dense shell. In this case, the X-ray 
emitting gas is located between the reverse shock and the 
swept-up outer shell, and is heated by the shocked ejecta 
via thermal conduction. Most of the mass in the hot region 
must come from mass-loading of denser interstellar clouds 
during the bubble expansion, and from the evaporation of 
part of the swept-up shell of interstellar medium, and its 
mixing with the lower-density, hotter jet material. 



The physical size of the X-ray hot spots is < 20 pc 
in radius, and the projected size of the whole bubble is 
~ 300 x 150 pc. From the estimated emission measures 
(Tables 3-5), we infer a mass of X-ray emitting gas ~ a 
few lOOM© in the hot spots (see also PSM10) and ~ a few 
lOOOM© in the cocoon (assuming a filling factor ~ 1). These 
values are several orders of magnitude higher than the mass 
that could have been carried out by the BH jet and winds 
over the source lifetime. But the mass of X-ray emitting gas 
is an order of magnitude less than the total mass of the 
swept-up interstellar medium; that is, the swept-up shell is 
not significantly depleted by evaporation into the hot re- 
gion, in agreement with the self-similar approximation of 
IWeaver et al.1 (1 19771 ). From the estimated mass and fitted 
X-ray temperatures, we conclude that the X-ray-emitting 
gas contains a thermal energy ~ 10 51 erg (hot spots) and 
~ 10 52 erg (whole cocoon). And we showed earlier that the 
total thermal energy ~ 10 53 erg and the energy carried by 
the synchrotron-emitting relativistic electrons ~ a few 10 50 
erg. The cooling timescale of the X-ray emitting gas in the 
cocoon is ~ 10 r yr, 100 times longer than the age of the 
source. This is consistent with an emitted X-ray luminosity 
~ 10 37 erg s" 1 even though power may have been trans- 
ferred from the jet to the ~ 0.3-1 keV component of the gas 
at an average rate ~ 10 39 erg s" 1 . By analogy with other 
shock-heated bubbles, we suggest that there may be even 
hotter but much less dense gas components, especially near 
the hot spots, whose hard X-ray emission emission would be 
too faint to be detectable in the 50-ks Chandra observation. 
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6 CONCLUSIONS 

We have carried out a multiband study of the shock-ionized 
bubble S26 in NGC 7793, which looks like a long-sought ana- 
logue of the Galactic jet source SS433/W50 but on an even 
grander scale (projected size « 300 x 150 pc). We showed 
that its structure is a scaled-down version of powerful FRII 
radio galaxies, with a core, radio lobes, X-ray hot spots and 
cocoon. It is the first time that all these elements have been 
found in a non-nuclear BH. We showed that the radio and 
X-ray hot spots are not spatially coincident: the X-ray hot 
spots are ~ 20 pc further out than the peak of the radio 
intensity in the lobes. This suggests that X-ray and radio 
emission come from different populations of radiating par- 
ticles. Based on our Chandra spectral analysis, we argued 
that the X-ray emission from the hot spots is most likely 
thermal. From the ATCA data, we showed that the radio 
emission from the lobes has a steep spectrum, consistent 
with optically-thin synchrotron emission. Over the rest of co- 
coon, the radio spectrum is flatter, suggesting an additional 
contribution from free-free emission; this is consistent with 
what we would expect from the measured H/3 line emission. 
A point-like radio core is not detected, but the radio spec- 
trum is flat or inverted in the proximity of the X-ray /optical 
core position; this may be interpreted as a more recent ejec- 
tion. However, deeper ATCA observations are needed (and 
scheduled) to test this suggestion. 

The total particle energy in the bubble is ~ 10 53 erg. 
Based on the measured radio flux and size of the bubble, and 
using standard equipartition relations for microquasar lobes, 
we estimated that the energy carried by the synchrotron- 
emitting relativistic electrons is a few 100 times less than 
the energy stored in protons, nuclei and non-relativistic elec- 
trons; non-relativistic particles provide most of the pressure 
to inflate the bubble. This system can give us important 
clues on how BHs at near-Eddington accretion rates trans- 
fer energy to the surrounding medium. The size and total 
energy content of the bubble are comparable to those found 
in some ULXs. However, here the core appears to be cur- 
rently X-ray faint (and was so also during the Einstein and 
ROSAT observations), while the jet is carrying a long-term- 
average power ~ 10 40 erg s- 1 (PSM10). We do not have 
any information on the long-term-average X-ray luminosity 
of the core, so we cannot exclude that it is similar to the me- 
chanical power. If the BH in S26 is of stellar origin, its super- 
Eddington jet power may force us to rethink the "canonical" 
scheme of BH accretion states. In Galactic BH transients, a 
collimated jet is present at accretion rates < a few percent 
of the Eddington rate (low/hard state). At higher accretion 
rates (~ 0.05-0.5 Eddington), the accretion flow usually 
collapses to a geometrically-thin, radiatively-efficient ther- 
mal disk, and the jet is quenched. At even higher accretion 
rates (above Eddington), high X-ray luminosity and power- 
ful mass-loaded outflows may coexist, but it is not known 
whether there can also be steady, collimated jets, and what 
their power is compared with the radiative power. S26 sug- 
gests that there can be collimated jets, and they may even 
dominate over the radiative output. The same scenario has 
been su ggested for some powerful FR II radio galaxies and 
quasars |Punslvll2007l ; llto et alj|2008j ). 
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APPENDIX A: SYNCHROTRON EMISSION IN 
THE MINIMUM-ENERGY CONDITION 

To estimate the minimum energy associated with the 
synchrotron-emitting cocoon, we assu me an energy range 
(7mm , 7max) for the relativistic electrons (|Poh l 1993; B icknelll 
2005), rather than a frequency range. Typical e mpirical val- 
ues of 7 m in ~ 1-10 (|Blundell fc Rawlingsll200oT ) and 7 max ~ 
10 4 -10 5 . For a steep spectrum, the minimum energy depends 
only very weakly on the high-energy cut-off. We introduce 
the following quantities: e e is the energy density in relativis- 
tic electrons; e v = (1 + k)e e is the energy density in rela- 
tivistic particles (electrons and protons), where fc is a free 
parameter; ea is the energy density in the magnetic field; 
etot = e p + en includes the energy in relativistic particles 
and magnetic field; e' tot is the total energy density including 
relativistic and non relativistic particles and the field. We in- 
troduce another free parameter n = etot/ttot to express the 
relative fraction of total energy stored in relativistic par- 
ticles (protons and electrons) plus field. The main reason 
why we distinguish between e to t and e' tnt is that there is 
solid observational evidence dCavagnolo et al1l2010l ; IPunslvl 
120071 ; iLeahv fc Gizanill200ll ; IWillott et alj|l999h that most 
of the energy in the lobes and cavities of radio galaxies is 
in low-energy electrons and other non-relativistic particles 
(i.e. r\ <C 1), and that the energy density of the magnetic 
field may be ~ 10-100 times less than the total energy den- 
sity e[ ot . Thus, using e' tot to derive a minimum-energy or 
equipartition criterion generally leads to very inaccurate es- 
timates for the synchrotron emission of a lobe. Instead, here 
we apply those criteria only to e to t, that is we assume that 
magnetic energy density is of the same order of magnitude 
as the relativistic particle energy density but much less than 
the total particle energy density. And within the relativis- 
tic energy density component, we use the parameter k to 
express the relative contribution of nuclei and electrons. 

Applying the minimum- energy condition leads, after 
some algebra |Bicknelll [20051 ). to this expression for the 
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minimum-energy magnetic field 
2 r p+l 



B 



ft(p,7min,7n 



C 

m, 

4/ (p+5) 



(1 + 

4/(p+5) / 3 \4/(p+5) 



4/(p+5) 



I e / \m e / V47T 



p+1 



(1 + fc) C~ x h 
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4/(p+5) 



5, 



4/(p+5) , 2(p-l)/(p+5) 



d 8/(p+5) r -12/(p+5) 
(Al) 



where the energy spectrum of the electrons is N(E)dE ~ 
E~ p dE, m e and e are the electron mass and charge, c the 
speed of light, I v the specific surface brightness S v the spe- 
cific flux at the observer's position, integrated over the whole 
cocoon, r c is the cocoon radius, and d is the distance to the 
source. We have assumed a filling factor of 1, for simplicity. 
The functions 

1 



h(p, 7min,7max) = 

C(P) 



p-2 

gp/2 



7 ( 2 .- p > 
' mm 
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(A2) 
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r(§ + if)r(f-^)r(f + i) 
r(f + I) 
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and F(z) is the Gamma function. The corresponding total 
(minimum) energy density is: 



etot.min = [(1 + k)e e + € B ] 



>+l 



B- 



2/io 



(A4) 



We now need to relate the energy density etot.min to 
the input jet power and size of the bubble. An approximate 
expression we could use is that the total energy (relativis- 
tic, non relativistic and field) is simply Vt. However, this is 
not entirely correct, because part of the injected energy is 
spent to inflate the bubble. From the self-similar solution of 
IWeaver et al.l |l977), we obtain a more accurate expression 
for the energy still available: 

3 5 



£tot = I^iT vtr ' 



and according to our definition of r], 
£tot S ilil 7 ^" 3 - 

In the minimum-energy approximation, from Eq. A4: 



3 5 



2^o I 4tt ll 77 
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P + 5 
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(A5) 
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(A7) 



and this value can now be substituted into Eq. Al. Fi- 
nally, the cocoon radius r c is obtained from the lWeaver et all 
1 19771 ) 's set of self-similar solutions (assuming a thin outer 
shell): 



r ~( 125 V /5 



1/5 



0.76 x 



Vt 3 
pa 



1/5 



(A8) 



and this expression is also substituted into Eq. Al. (Note 
that here it is the total jet energy Vt that determines the 
size of the bubble). 



Rearranging Eq. Al with such substitutions, we obtain: 

(p+5)/4 



S„ 



1.84 0.40 7? 
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p + 5 



„ 2\(P+5)/ 4 

2fj, Q e \ m 
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(A9) 



where the numerical values of h and C come from 
Eqs. (A2,3). For a spectral index a ~ —0.5 (p ~ 2), 
7min < 10 and 7 m ax ~ 10 5 , we have, in physical units: 

S„ » 640 (1+fc)" 1 77 7/4 P 39 13 t° 5 A n? 45 dT 2 us ' 5 mJy, (A10) 

where V39 is the jet power in units of 10 39 erg s _1 , t$ is the 
source age in units of 10 5 yr, ri\ is the interstellar number 
density in cm -3 , d is the source distance in Mpc, ^5 the 
observed frequency in units of 5 GHz; the numerical coef- 
ficient is not very sensitive to the choice of 7mm- We can 
obtain an analogous estimate for a spectral index a ~ —0.7 
(corresponding to p m 2.4), which is more often the case in 
radio lobes. In that case, 



S„ » 82 (1+fc) 



-1 1.85 ,t, 1.34,0.32 0.51 j-S — 0.7 t 

V V39 t 5 Hi d 1 v 5 mjy, (All) 



where we have fixed this time 7 m i n = 1. 

We can now compare these specific fluxes with the 
observations: S26 has a 5.5-GHz flux ~ 2 mjy, for a jet 
power ~ a few 10 40 erg s _1 , at a distance of 3.9 Mpc. This 



tells us that (1 + fc) 1 x 



V 



10" 



that is the energy 



stored in synchrotron-emitting relativistic electrons (a frac- 
tion rj/(l + k) of the total) is much less than the energy 
stored in relativistic protons and in non-relativistic parti- 
cles. We cannot determine the individual values of r\ and fc 
from this set of radio observations alone: only their combi- 
nation. If we use cosmic rays as an analogy (since they may 
be accelerated in jet and supernova shocks), we would ex- 
pect fc ~ 100. For plausible values fc ~ 10-100, 77/(1 + fc) ~ 
few 10~ 3 . 
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